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1.  Introduction 


The  demand  for  high  temperature  capacitor  dielectrics  in  power  electronics  is  increasing  with 

electrification  of  military  propulsion  and  weapons  systems.  Moving  toward  wide  bandgap 

semiconductors  (e.g.,  silicon  carbide  [SiC])  will  enable  operating  temperatures  above  140  °C, 

which  eases  thermal  management.  However,  such  systems  cannot  be  designed  efficiently  unless 

capacitors  are  available  that  can  operate  at  similarly  high  temperatures.  Metalized  polymer  film 

capacitors  have  the  advantage  of  self-healing,  which  allows  graceful  failure,  i.e.,  a  gradual  loss  in 

capacitance,  rather  than  catastrophic  failure  as  in  ceramic  or  film/foil  counterparts.  As  a  result, 

the  capacitor  dielectric  can  be  operated  at  an  electric  field  near  dielectric  breakdown,  thus 

achieving  a  higher  energy  density.  Presently,  the  widely  used  dielectrics  are  biaxially  oriented 

polypropylene  (PP)  and  poly(ethylene  terephthalate)  (PET),  which  are  low  cost,  available  in  thin 

film  (>0.5  pm  in  the  case  of  PET  and  >2.5  pm  for  PP),  and  have  a  relatively  high  breakdown 

field  (Ebd  >  550  MV/m).  Another  important  attribute  that  PP  offers  is  a  low  dissipation  factor 

(tan  5  --lx  10^)  that  is  independent  of  frequency.  In  comparison,  PET  is  rather  lossy 
—2 

(tan  5-1x10  ),  even  though  its  dielectric  constant  is  about  3.2,  which  is  50%  higher  than  PP. 

The  major  disadvantage  of  both  of  these  dielectrics  is  that  the  maximum  operating  temperature  is 
100  °C  for  PP  and  125  °C  for  PET.  Based  on  these  materials,  the  capacitors  occupy  about  25% 
of  volume  in  the  entire  system.  Therefore,  alternative  dielectrics  that  can  operate  at  high 
temperatures  are  needed  to  alleviate  thermal  management  and  improve  the  overall  energy  density 
(volume  efficiency). 

Some  of  the  high  temperature  polymers  that  are  available  in  commercial  thin  films  (10  pm  or 
less)  include  poly(phenylene  sulfide)  (PPS),  poly( ether-ether-ketone)  (PEEK),  and  poly( ether 
imide)  (PEI),  the  chemical  structures  and  thermal  properties  of  which  are  shown  in  Fig.  1 .  PP 
and  PET  are  also  included  for  comparison. 
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Among  the  three  high  temperature  polymers,  PEEK  has  potential  due  to  its  strong 
meehanieal/thermal  stability,  as  indicated  by  its  relatively  high  Tg  and  Tm,  as  shown  in  Table  1. 
Based  on  previous  results  from  electrical/dielectric  characterization,  PEEK  has  other  positive 
attributes.  One  of  them  is  a  relatively  low  dissipation  at  10  kHz  (the  typical  operating  frequency 
for  power  conditioning  type  applications)  for  temperatures  above  1 10  °C.  A  comparison  of  the 
loss  tangent  of  the  polymers  is  shown  in  Eig.  2.  Another  positive  attribute  that  PEEK  has  is  the 
rather  high  volumetric  resistivity  at  elevated  temperatures  for  fields  below  100  MV/m,  as  shown 
in  Eig.  3.  However,  the  volumetric  resistivity  decreases  significantly  as  the  field  is  increased 
further.  Another  limitation  of  PEEK  is  the  relatively  low  breakdown  field  (-400  MV/m) 
compared  to  the  other  polymers,  as  indicated  in  Eig.  4. 


Table  1  Thermal  properties  of  various  heat  resistant  polymers 


Polymer 

rp  a 

ro 

rp  b 

^  m 

(°C) 

Percent 

Crystallinity 

Dielectric 

Constant 

PPS 

120 

285 

~40 

~3 

PEEK 

150 

340 

~30 

3.1 

PEI 

220 

none 

none 

3.2 

PET 

75-80 

255 

~50 

~3.2 

PP 

<25 

170 

>60 

2.25 

“  Tg  denotes  glass  transition  temperature. 
^  Tjj,  denotes  crystal  melting  temperature. 


Fig.  2  Loss  tangent  at  10  kHz  as  a  function  of  temperature  for  various  polymers 
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Fig.  3  Volumetric  resistivity  as  a  function  of  electric  field  at  elevated  temperature  for  various  polymers 
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In  order  to  gain  a  fundamental  understanding  of  the  underlying  eauses  and  mechanisms  that  lead 
to  high  field  conduction  and  low  breakdown  field,  we  aim  to  study  the  electronic  structure  of 
PEEK  and  investigate  the  effects  of  chemical  impurities,  such  as  water  molecules,  and  defects, 
such  as  chain  ends,  on  the  band  gap  through  ah  initio  calculations.  In  particular,  our  goal  is  to 
determine  if  any  intermediate  energy  states  (impurity  states)  are  introduced  in  the  band  gap.  In 
this  report,  progress  on  the  study  of  electronic  structure  and  band-gap  calculation  of  pristine 
PEEK  is  summarized. 


2.  Computational  Methodology 


All  quantum  mechanical  calculations  were  performed  using  density  functional  theory  (DET) 

1  2 

methods  with  a  mixed  Gaussian  and  plane -wave  (GPW)  ’  approach  as  implemented  in  the 
CP2K  package.  The  exchange-correlation  functional  was  described  by  the  generalized  gradient 
approximation  (GGA)  with  the  spin-polarized  functional  developed  by  Perdew-Burke-Emzerhof 
(PBE)."^  Grimme’s  D3  pair-wise  dispersion  correction  (DFTD3)^  was  employed  to  account  for 
van  der  Waals  interaction.  Wavefunctions  were  expanded  in  valence  triple-zeta  plus  polarization 
(TZVP)  Gaussian  basis  sets^  with  an  auxiliary  plane-wave  basis  with  a  cutoff  energy  of 
900  Rydberg  (Ry).  The  plane-wave  cutoff  of  900  Ry  was  determined  by  converging  the  energy 
with  respect  to  the  plane-wave  cutoff  to  1  meV  per  atom.  Core  electrons  were  modeled  by  the 
Goedecker-Teter-Hutter  (GTH)^  scalar  relativistic  norm-conserving  pseudopotentials  with  1,  4, 
and  6  valence  electrons  for  hydrogen  (H),  carbon  (C),  and  oxygen  (O),  respectively.  Brillouin 
zone  integration  was  performed  with  a  reciprocal  space  mesh  consisting  of  only  the  E-point. 

The  initial  geometry  for  all  calculations  was  constructed  based  on  the  crystal  unit  cell  parameters 

o 

of  the  6-ring  (i.e.,  2  repeating  units  of  PEEK)  model  developed  by  lannelli.  In  order  to  reduce 
errors  due  to  periodic  approximations,  the  crystal  unit  cell,  which  comprised  of  2  chains  of  2 
repeating  units,  was  doubled  along  the  y-axis  to  create  a  supercell  of  4  chains,  totaling 
272  atoms.  The  supercell  is  repeated  using  periodic  boundary  conditions  to  simulate  an  infinite 
system.  Each  supercell  was  optimized  by  allowing  the  lattice  vectors  and  lattice  angles  to  shift. 
Initial  geometry  optimizations  were  performed  at  the  PBE-I-DFTD3/TZVP  level  of  theory  with  a 
convergence  criterion  of  0.02  eV/A  or  less  for  the  force  on  each  atom. 

After  geometry  optimization,  the  ground-state  (i.e.,  temperature  at  0  K)  electronic  structure  of 
PEEK  was  determined,  from  which  the  electron  density,  density  of  states  (DOS),  projected 
density  of  states  (PDOS),  and  band  gap  can  be  determined.  Initial  calculations  using  the  PBE- 
GGA  exchange-correlation  functional  failed  to  determine  the  band  gap  accurately;  therefore,  we 
performed  single-point  calculations  using  the  hybrid  functional  developed  by  Heyd-Scuseria- 
Emzerhof  (HSE06)^  on  the  optimized  structure  described  previously.  This  approach  was 
validated  on  the  band  gap  of  polyethylene  with  1 00%  agreement. 
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For  large  systems  (e.g.,  PEEK)  deseribed  with  high-quality  basis  sets,  ealeulation  of  the  Hartree- 
Eoek  exaet-exehange  (HEX)  eomponent  in  the  exehange-eorrelation  funetional,  sueh  as  that 
ineorporated  in  the  HSE06  functional,  is  computationally  demanding.  To  speed  up  the 
calculations,  we  employed  the  auxiliary  density  matrix  method  (ADMM),^°  which  has  been 
shown  to  outperform  the  standard  HEX  implementation  by  a  factor  of  20  with  excellent 
performance  and  good  accuracy.  In  this  work,  the  ADMM  was  employed  with  the  TZVP  GTH 
primary  basis  set  and  PITS  auxiliary  fit  basis  set. 

To  reassure  that  the  computed  structure  represents  the  system  under  investigation,  we  compared 
the  computed  x-ray  diffraction  (XRD)  pattern  with  our  measurement.  The  XRD  calculations 
reported  in  this  work  employ  the  method  implemented  by  Krack  and  co-workers^  ^  in  the  CP2K 
code,  in  which  the  elastic  or  coherent  x-ray  scattering  intensities  I(q)  within  the  first  Bom 
approximation  for  systems  where  no  angular  resolution  exists  are  given  by 

Kq)  =  {\F(.q)n  (1) 

where 

F(q)  =  J  n(r)e‘‘i*‘dr,  (2) 

is  the  Pourier  transform  of  the  1 -electron  density  n(r),  r  is  the  position,  q  is  the  scattering  vector, 
and  the  average  in  I(q)  is  over  all  q  vectors  with  the  same  modulus  q  in  the  statistical  ensemble. 
F(q)  was  calculated  using  the  all-electron  Gaussian  augmented  plane-wave  (GAPW)  ’  method, 
in  which  the  Kohn-Sham  orbitals  are  expanded  in  plane-waves  with  the  electronic  charge 
density  divided  into  a  smooth  part  that  is  expanded  in  plane-waves  and  an  atom-centered,  more 
rapidly  varying  part,  which  is  expanded  in  Gaussian  functions.  As  pointed  out  by  Krack,  the  all¬ 
electron  method  was  employed  instead  of  the  plane-wave-pseudopotential  approach  since  the 
latter  tends  to  give  an  incorrect  description  of  the  electronic  charge  density  close  to  the  nuclei. 
The  Becke,  Eee,  Yang,  and  Parr  (BEYP)  ’  semilocal  form  of  the  exchange-correlation 
functional  was  used  in  all  XRD  calculations.  Other  functionals  such  as  PBE,  B3EYP^"^,  and 
HSE06  were  also  compared,  and  no  significant  difference  from  BEYP  was  found. 


3.  Results 


Three  initial  geometries  of  PEEK  were  constmcted  to  determine  the  best  representation  for  our 
system.  As  shown  in  Pig.  5,  Geometry  I  was  created  by  adapting  the  unit  cell  report  by  lannelli 
without  any  manipulation  of  the  chains,  while  Geometry  II  was  constmcted  by  shifting  the 
second  chain  in  the  unit  cell  by  half  of  the  unit  cell  length  relative  to  first  one  along  the  z-axis 
and  Geometry  III  was  generated  by  rotating  the  second  chain  90°  around  the  z-axis  relative  to  the 
first  one.  All  three  geometries  contain  two  unit  cells  along  the  y-axis.  Table  2  summarizes  the 
supercell  parameters  of  the  three  optimized  geometries  and  the  deviation  denoted  by  A  from  the 
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literature  values.  ’  As  indieated  in  Table  2,  the  optimized  Geometry  III  gives  the  best  agreement 
with  the  literature  values.  The  finding  is  further  supported  by  the  computed  XRD  patterns,  which 
are  shown  in  Fig.  6.  These  results  serve  as  an  assurance  that  the  ground-state  (i.e.,  temperature  at 
0  K)  Geometry  III  computed  by  DFT,  which  agrees  well  with  experimental  measurements,  may 
be  used  for  future  calculations. 


Fig.  5  Illustration  of  the  3  optimized  supercells.  Chains  in  a  lighter  color  represent  ones  in  the 
background. 

Table  2  Comparison  of  the  structure  parameters  of  the  3  optimized  supercells  with  literature 


Literature 

I 

II 

III 

a  /  A  (A) 

7.78  /- 

9.30/  +7.52 

1.11 1-0.01 

7.58/ -0.2 

b  /  A  (A) 

11.84  /- 

11.74 /-tt7 

12.46/  +0.62 

11.89/  +0.05 

c  /  A  (A) 

29.92  /- 

30.21  /  +0.29 

29.83  / -0.09 

30.05/  +0.75 

a  /  A  (deg.) 

90.0  /- 

91.16/  +7.7(5 

89.96/ -0.04 

89.89/ -0.77 

(5  /  A  (  deg.) 

90.0  /- 

118.43/  +2&  45 

89.80/ -0.2 

90.04  /  +0.04 

Y  /  A  (deg.) 

90.0  /- 

80. 19/ -9.  <57 

90.10/  +0.7 

90.02  /  +0.02 

Volume  /  A  (A^) 

2756 /- 

2851  /  +95 

2888  /  +752 

2708  /  -48 
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Fig.  6  Comparison  of  the  computed  XRD  patterns  for  the  3  optimized  geometries  with  the 
experiment 

Geometry  III  was  used  as  the  basis  for  ealeulating  the  DOS  of  PEEK,  which  is  shown  in  Eig.  7a, 
from  which  the  band  gap  was  determined  as  4.21  eV,  which  is  33%  higher  than  the  experimental 
value  of  3. 1  eV  calculated  based  on  the  ultraviolet  (UV)  cut-off  wavelength  of  390  nm  shown  in 
Eig.  8  by 
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Eg  =  hdX, 


(3) 


where  Eg  is  the  band  gap,  h  the  Planek  eonstant  in  eV-s,  c  the  speed  of  light  in  nm/s,  and  A.  the 
UV  cutoff  wavelength  in  nm.  The  discrepancy  is  the  consequence  of  the  approximation  within 
DFT  for  describing  the  exchange-correlation  energy  of  electrons,  for  which  the  HSE06  hybrid 
functional  was  used  in  this  work.  To  obtain  better  agreement  for  the  band  gap  value  for  future 
calculations,  a  systemic  study  was  conducted  for  the  (a,  w)  parameter  space  of  the  HSE 
functional,  where  a  is  the  exchange  mixing  coefficient  that  determines  the  amount  of  exact 
exchange,  and  ro  is  the  screening  parameter  controlling  short-range  interaction.  The  band  gap 
converged  to  around  experimental  value  at  a  =  0.093  and  ro  =  0.1 1  bohr“*.  The  corresponding 
DOS  is  shown  in  Eig.  7b;  whereas,  the  one  shown  in  Eig.  7a  was  computed  from  a  =  0.25  and 
CO  =  0.1 1  (these  parameters  are  based  on  the  version  of  Geometry  III  that  was  forced  to  be 
orthogonal).  Euture  HSE  calculations  will  be  performed  with  the  updated  parameters. 


-25  -20  -15  -10  -5  0  5  10  15  20 

Energy  (eV) 


Fig.  7  Density  of  States  of  PEEK  based  on  the  optimized  Geometry  III: 
a)  from  standard  HSE06  hybrid  functional  with  a  =  0.25  and  co  = 
0.1 1  bohr  '  and  b)  from  HSE  with  a  =  0.093  and  co  of  0.1 1  bohr  ' 
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Fig.  8  UV  absorption  spectra  of  amorphous  and  semi-crystalline  PEEK.  PP  and  PEI  are  included  for 
comparison 


4.  Conclusion 


The  correct  geometry  of  the  unit  cell  and  a  methodology  for  computing  accurately  the  band  gap 
of  PEEK  were  developed.  Dispersion-corrected  DFT  was  used  to  optimize  various  starting 
geometries,  and  the  correct  one  was  determined  based  on  a  comparison  with  experimental 
crystallographic  data.  The  band  gap  was  determined  from  the  PDOS,  which  is  calculated  using 
the  HSE06  hybrid  exchange-correlation  functional.  The  original  formulation  of  the  HSE06 
functional  overestimates  the  band  gap  of  PEEK  by  33%  in  comparison  to  the  experimental  value. 
Therefore,  a  systematic  study  of  the  parameter  space  of  the  functional  was  conducted  by  varying 
the  exchange-correlation  mixing  and  screening  coefficients  to  obtain  a  functional  form  that  best 
reproduces  the  experimental  value  of  the  band  gap. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


ADMM 

auxiliary  density  matrix  method 

BLYP 

Becke,  Eee,  Yang,  and  Parr 

C 

carbon 

DFT 

density  functional  theory 

DFTD3 

D3  pair-wise  dispersion  correction 

DOS 

density  of  states 

GAPW 

Gaussian  augmented  plane-wave 

GGA 

generalized  gradient  approximation 

GPW 

Gaussian  and  plane-wave 

GTH 

Goedecker-T  eter-Hutter 

H 

hydrogen 

HFX 

Hartree-Eock  exact-exchange 

HSE 

Heyd-Scuseria-Emzerhof 

0 

oxygen 

PBE 

Perdew-Burke-Ernzerhof 

PDOS 

projected  density  of  states 

PEEK 

poly(  ether-ether-ketone) 

PEI 

poly( ether  imide) 

PET 

poly( ethylene  terephthalate) 

PP 

polypropylene 

PPS 

poly(phenylene  sulfide) 

Ry 

Rydberg 

SiC 

silicon  carbide 

TZVP 

variance  triple-zeta  plus  polarization 
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uv 

XRD 


ultraviolet 
x-ray  diffraction 
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